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Abstract. The paper is centered around a new proof of the infini- 
tesimal rigidity of convex polyhedra. The proof is based on studying 
derivatives of the discrete Hilbert-Einstein functional on the space of 
"warped polyhedra" with a fixed metric on the boundary. This ap- 
proach is in a sense dual to using derivatives of the volume in order to 
prove the Gauss infinitesimal rigidity of convex polyhedra, which deals 
with deformations that preserve face normals and face areas. 

In the spherical and in the hyperbolic-de Sitter space, there is a 
perfect duality between the Hilbert-Einstein functional and the volume, 
as well as between both kinds of rigidity. 

We also discuss directions for future research, including elementary 
proofs of the infinitesimal rigidity of hyperbolic (cone-)manifolds and 
development of a discrete Bochner technique. 
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1. Introduction 



1.1. Infinitesimal rigidity of convex polyhedra. Among various con- 
cepts of rigidity, the one that we deal with here is the infinitesimal rigidity. 
Let P C be a compact convex polyhedron; assume for simplicity that all 
faces of P are triangles. Then an infinitesimal deformation of P is an as- 
signment of a vector qi to each vertex pi. Let P{t) be a family of polyhedra 
with the same combinatorics as P and vertices Pi{t) = Pi + tqi. An infini- 
tesimal deformation is called isometric, if the edge lengths of P{t) remain 
constant in the first order of t at t = 0. There always exist trivial isomet- 
ric infinitesimal deformations that are restrictions of infinitesimal isometrics 



Every convex polyhedron P is infinitesimally rigid, that is 
every isometric infinitesimal deformation of P is trivial. 

This theorem was first stated by Dehn in [15]. Legendre-Cauchy's argument 
from [To] easily carries over to the infinitesimal setting. Other proofs are 
given in [i m Ea SSI Ea EH]. 

Yet another proof is presented in this paper. In some aspects it is similar 
to Schlenker's proof [l6] and to the Pogorelov-Volkov proof of global rigidity 
[39l [53] . The actual goal of this paper is to describe the general framework 
into which this argument fits and to indicate possible further developments. 

1.2. The approach. Instead of deforming the embedding P C M'^ by mov- 
ing the vertices, we deform the metric in the interior of P. We choose a 
point pq G int P and subdivide P into pyramids with po as the apex and 
faces of P as bases. By denoting = \\po — Pi\\, we start to vary lengths rj, 
while leaving the lengths of boundary edges constant. As a result, the total 
dihedral angles Ui around interior edges popi may become different from 27r, 
so that the polyhedron P becomes what we call a warped polyhedron. 

Can one vary the lengths rj so that the angles Ui remain constant in the 
first order? One can always do this by moving the point ro inside P. It 
turns out that if these are the only possibilities, then the polyhedron P is 
infinitesimally rigid: infinitesimal rigidity of P is equivalent to 



where Kj = 27r — Wj is the curvature of the edge poPi- Our proof of the 
infinitesimal rigidity of convex polyhedra goes by determining the rank of 
the above Jacobian. 

An immediate generalization of this approach is to consider an arbitrary 

triangulation of the polyhedron P and the Jacobi matrix (^^^^ > where rij 

and the length, respectively, curvature of the edge PiPj, and among 

the points pi some are vertices of P, some lie on its edges or faces, and 
some in the interior. (As before, lengths of boundary edges are assumed 



of ]R3. 



(1.1) 
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constant.) Then infinitesimal rigidity of P is equivalent to 



where m is the number of interior vertices, and n is the number of vertices 
interior to faces of P. In our joint paper with Jean-Marc Schlenker |28] . we 
have shown that the above Jacobian has exactly m positive eigenvalues. In 
particular, the matrix has negative spectrum if all points pi lie at vertices 
or on the edges of P. 

1.3. Hyperbolic manifolds with boundary and hyperbolic cone- 
manifolds. Similar to the previous subsection, consider a triangulation of 

a compact closed hyperbolic 3-manifold and the Jacobi matrix of 

edge curvatures with respect to edge lengths. 



Conjecture 1. The matrix y'Q:p^j has corank 3n and exactly n positive 
eigenvalues, where n is the number of vertices in the triangulation. 

The first part, corank equals 3n, is known. It follows from the infinitesimal 
rigidity of compact closed hyperbolic manifolds in dimension greater than 2, 
the so-called Calabi-Weil rigidity, [8l[57]. Therefore, if the arguments used 
in this paper can be extended to prove Conjecture [H this would yield an 
elementary proof of Calabi-Weil rigidity. Besides, such extension should lead 
to a proof of the following conjecture, communicated to me by Jean-Marc 
Schlenker who proved a smooth analog in j45j . 

Conjecture 2. Compact hyperbolic manifolds with convex polyhedral bound- 
ary are infinite simally rigid. 

More generally, consider a compact closed manifold glued from hyperbolic 
simplices so that the total dihedral angles around the edges may be different 
from 271. It is called a (triangulated) hyperbolic cone-manifold. Due to 
works of Hodgson-Kerckhoff [23], Mazzeo-Montcouquiol [32], and Hartmut 
Weiss [58], it is known that compact closed hyperbolic manifolds with cone 
angles less than 27r are infinitesimally rigid, that is cannot be deformed so 
that their cone angles remain constant in the first order. Again, this can be 
reformulated in terms of the rank of the Jacobian of the map r i— )> k, and a 
generalization of our method would yield an alternative, elementary proof. 

Note that the condition on cone angles is similar to the convexity condition 
for polyhedra. There exist infinitesimally flexible non-convex polyhedra as 
well as infinitesimally flexible cone-manifolds with some cone angles greater 
than 27r, see [HIETj. 

A study of deformations of Euclidean cone-manifolds from a similar point 
of view is undertaken in |21] . 
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1.4. Discrete Hilbert-Einstein functional. Here we sketch our proof of 
dimker = 3 for a star-like triangulation of a convex polyhedron P. 
Consider the function 

HE(r) = ^ TiKi + ^ ^ijXij, 

where ^ij is the length of, and Xij is the exterior dihedral angle at an edge 
PiPj of P. We call it the discrete Hilbert-Einstein functional, as it is the 
discrete analog of twice the total scalar curvature of P plus half of the total 
mean curvature of dP. The Schlafli formula implies that = Kj, therefore 

dKi _ ^2 HE 

dvj dvidrj ' 

in particular the left hand side is symmetric in i and j. This allows to 
reformulate the infinitesimal rigidity in the following way: 

(1.2) If HE" = 0, then r is a trivial variation of r. 

Here HE" = ■ f^.™. 'f'i'f'j = f^'i'i'i is the second derivative of HE in the 
direction r, and a trivial variation can be defined as the one that preserves 
boundary dihedral angles: Xij = 0. 

In order to prove p.2|) . we derive formula (j3.12|) for HE". This formula 
implies that HE" is non-positive, and vanishes only if r is trivial. But as 
HE" vanishes by assumption, the variation r must be trivial. 

This argument is reminiscent of Koiso's proof of the infinitesimal 
rigidity of Einstein manifolds under certain assumptions on the curvature 
operator. Koiso uses integration by parts to obtain two formulas for the 
second derivative of the Hilbert-Einstein functional. This yields an equation 
with zero on one side, while on the other side one has a non-positive quantity 
that vanishes only if the deformation is trivial. Note that our formula (j3.12p 
is obtained from HE" = kiVi by kind of discrete integration by parts. 

Koiso's proof is an example of application of Bochner's technique: a sec- 
ond order differential operator is expressed as the sum of Laplacian and of 
a non-negative 0-th order operator. Thus our proof should be a particular 
manifestation of a discrete Bochner technique. 

1.5. Discrete Bochner technique. Einstein manifolds in dimension 3 are 
manifolds of constant sectional curvature, therefore Koiso's theorem contains 
Calabi-Weil rigidity of hyperbolic 3-manifolds as a special case. Weil also 
uses Bochner's technique; his approach is related to that of Koiso in the 
same way as moving vertices of a polyhedron P is related to deforming the 
metric inside P. The infinitesimal rigidity of hyperbolic cone-manifolds with 
cone angles less than 2-k (see Subsection II. 3|) is proved by extending Weil's 
arguments with the help of Cheeger's Hodge theory for singular spaces. 

It will be just natural if the infinitesimal rigidity of cone-manifolds can 
be reproved by developing a discrete Bochner technique for this situation. 
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Such a proof would not only be elementary, it would also provide a discrete- 
geometric counterpart to the original argument of Hodgson, Kerckhoff, and 
others. 

Another theorem where the need for a discrete Bochner technique is felt, is 
Cheeger's discrete analog [TT] of Bochner-Gallot-Meyer vanishing theorem: 

If in a Euclidean cone-manifold of dimension d all cone angles 
are less than 27r, then it is a real homology d-sphere. 

Again, the known proof uses Hodge theory for singular spaces. 

It should be mentioned that Forman developed a combinatorial Bochner 
technique |20| that takes into account the combinatorics of a simplicial com- 
plex (or, more generally, cell complex), but not its geometry, expressed by 
edge lengths. 

One of the components for a discrete Bochner technique should be the 
discrete Hodge theory based on Whitney forms, see e. g. [16]. The other 
component presumes some sort of discrete Riemannian geometry, still to be 
found. 

1.6. Volume derivatives and the Alexandrov-Fenchel inequality. If 

one goes into details of the proof of dimker = 3 sketched in Subsection 

11.41 then one sees that it resembles very much a known proof of 

(1.3) dimker (l^i) =3, 

where hi are lengths of perpendiculars dropped to faces of a convex polyhe- 
dron Q from the origin, Ai are face areas, and the polyhedron is deformed 
by varying hi while keeping the directions of face normals fixed. 

In fact, equation (jl.Sp . with d on the right hand side, also holds for 
convex d-dimensional polyhedra. This is a key lemma in the proof of the 
Alexandrov-Fenchel inequalities. Alexandrov-Fenchel inequalities also imply 
that the Jacobian in (jl.3p has exactly one positive eigenvalue (cf. the end 
of Subsection II. 2p and describe a part of the positive cone. 

It is geometrically clear that Ai = where Yol{h) is the volume of a 

polyhedron with given face normals and with support numbers hi. Therefore 
equation (jl.Sp computes the rank of the second derivative of the volume, thus 
completing the analogy with Subsection 11.41 

1.7. Duality between the volume and Hilbert-Einstein functional. 

The analogy described in the previous subsection culminates in a striking 
identity 

92 HE _ a^Vol 
dridvj dhidhj ' 

where the polyhedron Q on the right hand side is polar dual to the poly- 
hedron P on the left hand side. We were able to prove the above identity 
only by a direct computation. Note that functions HE(r) and Vol(/i) have 
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different nature: the former uses inverse trigonometric functions while the 
latter is polynomial. 

The situation is nicer in spherical and hyperbolic geometry. For a spher- 
ical convex polyhedron P and its polar dual P* we have by [33] 




This implies S{P) + S*{P*) = 27r^, where S and S* are functionals on 
the space of warped spherical polyhedra. Functionals S and S* are the true 
analogs of HE and Vol, respectively; in particular their variational properties 
are similar to those of their Euclidean counterparts. For more details, see 
Subsections 14.21 and 14.31 The hyperbolic case is similar, but the polar dual 
P* of a hyperbolic polyhedron P lives in the de Sitter space, see Subsection 



1.8. Gauss infinitesimal rigidity and Minkowski theorem. While equa- 
tion (jl.ip is equivalent to the infinitesimal rigidity of a polyhedron, equation 
(jl.3p implies what we call Gauss infinitesimal rigidity: 

If the support numbers hi of a convex polyhedron vary in 
such a way that the face areas remain constant in the first 
order, then the polyhedron undergoes a parallel translation. 

This theorem was stated and proved by Alexandrov ^ Chapter XI]. Minkow- 
ski proved what we would call global Gauss rigidity: two convex polyhedra 
with the same face normals and face areas differ by a parallel translations. 
Existence is asserted in the Minkowski theorem: 

Given unit vectors Ui that span and positive numbers Ci 
such that ^ - CiVi = 0, there exists a convex polyhedron with 
face areas Ci and outward face normals fj. 

More generally, Minkowski existence and uniqueness theorem holds in M'^. 

1.9. Alexandrov theorem on existence of a convex polyhedron with 
a given metric on the boundary. A counterpart to the Minkowski the- 
orem in dimension 3 is the Alexandrov theorem 0: 

Given a Euclidean cone-metric g on the 2-sphere with all 
cone angles less than 27r, there exists a convex polyhedron in 
with g as the intrinsic metric on the boundary. 

In a joint paper [6] with Alexander Bobenko, we gave a new proof of the 
Alexandrov theorem. The polyhedron P is obtained by constructing a family 
of warped polyhedra -Pt,0 < i < 1 where the curvatures Kj tend to as t 
tends to 1. The local existence of such family is based on the following 
property of the Jacobian of the map r ^ k: 




Yo\{P) + -Y^hjh, +Yo\{P*) 



= vr' 



.2 



dimker ( -- — 




0, if < < 5i for all i. 
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Here 6i is the angular defect of the i-th cone point on (S'^,g). For the 
global existence one has to make sure that polyhedra don't degenerate in 
the process of deformation. Note that the triangulation of the boundary of 
P may change in the process of deformation. See also Subsection [ 



1.10. Acknowledgements. I would like to thank Frangois Fillastre and 
Jean-Marc Schlenker for interesting discussions and useful remarks. 

2. Gauss rigidity of convex polyhedra 

2.1. The theorem. A set {Qt)o<t<e of compact convex polyhedra is called 
a linear family, if each Qt is obtained from Q = Qo hy parallelly translating 
the planes of the faces, with the translation vector of each face depending 
linearly on t. 

In this section, a proof of the following theorem is presented. 

Theorem 2.1. Assume that the area of each face in the linear family (Qt) 
is constant in the first order oft. Then all Qt are translates of Q. 

This theorem is proved by Alexandrov in jU Chapter XI] by two different 
methods. The proof given here is essentially the second proof of Alexandrov, 
but presented in a self-consistent way. Our purpose is to reveal that the proof 
is based on certain variational properties of the volume, in order to make 
the relationship with Section [3] more straightforward. 

Let us introduce some notations. Let Fi, F2, . . . , Fn be the faces of Q, 
and let z/j be the outer unit normal to the face F^. Then the equation of the 
plane of the face Fi is 

span(Fj) = {x G I (x, Ui) = 

where /i^ is the signed distance from the coordinate origin G to the 
plane spanned by Fj. The numbers (h^) are called support parameters of 
the polyhedron Q. We have 

Q = {x G I (x, Ui) < /jO, i = 1, . . . , n} =: Q{h°). 

We fix the directions of outer normals, and let the support parameters 
vary. A linear family of polyhedra is Q{h^ + tu). A polyhedron Q{h) with 
hi = h^ + (a, i^i) is a parallel translate of Q by a G M^. 

Let Ai be the area of the face F^. Denote by DAi{u) the derivative of Ai 
in the direction u: 

OA 
dhi 



DAM := I 



Ai{h + tu) = ^ ^T^Ui. 



t=o 

We will usually put h in place of u and then write 

DAi{h) =: Ai. 

Thus the statement of Theorem 12.11 can be rewritten as 

If h € is such that at h = h^ we have Ai = for all i, 
then there exists a G M'^ such that hi = (a, z/j) for all i. 
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It is easy to see that the space of "trivial deformations" hi = (a,fj) has 
dimension 3. Hence the theorem can also be reformulated as 
The Jacobian of the map (hi) (Ai) has corank 3; 



dim ker 



dhj 



h=hO 



2.2. The approach. Recall that we put 

Q{h) := {x € I (x, i/j) < hi, i = l,. . . , n}. 

The set Q{h) is a convex polyhedron for all h G M", but it may have less 
than n faces or even be empty. Let U C M" be a neighborhood of h^ such 
that Q{h) has n faces for all /i G C/. Consider the function 

Vol: U ^R, 

where Vol(/i) is the volume of the polyhedron Q{h). 

Lemma 2.2. The function Vol is continuously difjerentiable on U with 

where Ai{h) is the area of the i-th face of the polyhedron Q{h). 

Proof. The equation (|2.ip is geometrically obvious: as we shift the plane of 
the i-th face by e, we glue to (or cut from) Q a convex slice of thickness e. 
One side of the slice has area Ai, the other side has area Ai + 0(e). Hense 
the volume of the slice is eAi + o(e), and (j2.ip follows. 

Since its partial derivatives Ai are continuous, the function Vol is contin- 
uously differentiable. □ 

Equation (12.ip implies that the Jacobian of the map {hi) ^ [Ai) equals 
the matrix of the second differential of the function Vol. Here by the second 
differential we mean a symmetric bilinear form 



D'^Yo\{u,v) := 4- 
dt 



a^voi 



t=o ij dhidhj 



UiV 



This yields the following reformulation of Theorem 12.11 

The second differential of the volume at h = h^has corank 3: 

dim ker (D^^, Vol) = 3. 

In order to prove Theorem 12.11 we will compute the second variation of 
the volume in two different ways and compare the formulas obtained. Here 
the second variation Vol" is the quadratic form associated to Vol: 

Vol" := D^Yol{h,h). 

Equivalently, 



Yol{h + th). 
t=o 
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2.3. Orthoscheme decomposition. Recall that Fi is the face of Q{h) with 
outer normal fj. If Fi and Fj share an edge, denote this edge by Fij. 

Let qi be the foot of the perpendicular dropped from G to the plane 
spanned by Fi. For every pair of adjacent faces, drop perpendiculars from 
the points qi and qj to the line spanned by Fij. Their common foot will be 
denoted by qij = qji. Finally, denote by qijk the vertex of Q{h) where faces 
Fi, Fj, and Fk meet. 

Denote by hij the signed length of the segment qiqij, the sign being posi- 
tive if qi lies on the same side from the edge Fij as the polygon Fi . Similarly, 
let hijk be the signed length of the segment qijqijk- See Figure [H 



From now on, we assume that the polyhedron Q{h^) is simple, i. e. that 
at each of its vertices meet exactly three faces. Then, by choosing a neigh- 
borhood U of hP appropriately small, we can ensure that all polyhedra Q{h) 
with h are combinatorially isomorphic (under identification of faces with 
equal outward normals). The case of a non-simple Q{hP) is a bit subtle, and 
we explain in Subsection 12.61 how it can be treated. 

Under this assumption, the functions 



are linear. If Fi n Fj = or Fi r\ Fj r\ F^ = 0, then the corresponding 
functions hij and hijk ^ire not defined or can be put identically zero. 

Lemma 2.3. For all h ^ U , we have 




Figure 1. Definition of hij and hijk- 



hij : [/ ^ M, 
hijk'- U ^ M., 



(2.2) 




Proof. Denote by lij the length of the edge Fij- Then we have 
(2.3) (.ij = hijk ~\~ hiji. 
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where qijk and qiji are the endpoints of Fij. By substituting this into 

(2.4) ^^ = IY1 ^v^^^ 

j 

and substituting the result into 

(2.5) Vol(/i) = i^/i^Ai, 

i 

we obtain ([221) ■ 

Alternatively, the right hand side of (|2.2p can be seen as the sum of signed 
volumes of the orthoschemes OqiQijqijk- D 

Note that since hij and hij^ are linear functions of h, formula (12. 2p ex- 
presses Vol : [/ — >■ M as a third degree homogeneous polynomial in variables 
(hi). 

2.4. First and second variations of the volume and of the face areas. 

We start with an auxiliary lemma. 

Lemma 2.4. For an arbitrary variation h € M" we have 

(2.6) hihij + hjhji = hihij + hjhji, 

(2.7) hijhiji^ + hif^hii^j — hijhijj^ + ^ikhikj- 

Here hij and hjjk denote derivatives in the direction of h. 

Proof. Consider the quadrilateral OqiqijQj, see Figure O The angle ifij at 
equals the angle between the normals and z/j, and is therefore constant. 
On the other hand, we have 

hi- hi 

(2.8) (fij = arctan — ^ + arctan 

/ (/^ j 

By differentiating in the direction of /i, we obtain 

1 hij hi hij hi 1 hjihj hjihj 

hijhi hijhi hjihj hjihj 

Observing that hf + hfj = h'j + h'j- , we arrive at ()2.6p . 

If some of the lengths hi, hj, hij, hji are negative, the formula (j2.8p still 
remains valid. Care has to be taken when hi = or hj = 0. One can avoid 
these difficulties by noticing that both sides of (j2.6p are linear functions of 
hi and hj; so in order to prove ()2.6p for all h it suffices to check it for all h 
in some open set. By letting h vary over a small neighborhood of the point 
(1, 1, . . . , 1), we ensure that all segments on Figure [2] have positive length. 
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Figure 2. The quadrilateral Oqiqijqj- 

(It does not matter whether (1, 1, . . . , 1) € U, as the functions {hij) can be 
extended to R" by linearity.) 

Equation (12. 7p is proved in the same way. For the last part of the argu- 
ment, one should note that (hijk) are linear functions of (hij). 

Alternatively, the lemma can be proved by using explicit formulas ex- 
pressing hij in terms of hi and hj , see Figure [3l 




a — b cos 7 
sin 7 



Figure 3. Consider orthogonal projection on the side a. 

□ 

Lemma 2.5. The first variations of the face areas are given by 
(2.9) Ai = ^ ^ hijhijf^ = ^ ^ hijhijf^. 

j,k j,k 
Proof. Equations (j2.4p and ()2.3p imply 

Ai = o y^^jhijhjjk + hjjhjjk) 

~ 2 ^ ^ {hijhiji^ -\- hi^hii^j) -\- — ^ ^ (hijhij^ -\- hn^hn^j), 
{j,k} {j,k} 

where k} denotes the sum over unordered pairs of j and k. By using 
(|2.7p and converting back to the sum over ordered pairs, we obtain (12. Dh . □ 
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Lemma 2.6. The first variation of the volume is given by any of the fol- 
lowing formulas. 

(2.10) YoY =^hiAi 

i 

(2.11) vor = ^^M*. 

i 

Proof. From (|2.5p we have 

(2.12) Vol- = ^^M, + ij^Mi. 

i i 

By using (j2.4p . regrouping, then using ()2.6p . regrouping back, and finally 
applying Lemma 12.51 we obtain 

2 ^ ^ hiAi — ^ ^ /li ^ ^ hij£ij — ^ ^ hihijhijig — ^ ^ {hihij + hjhji)hijk 

i i j {i,j},k 

— ^ ^ (hihij -\- hjhji^h^ji^ — ^ ^ hihijhijk — ^ ^ ^ ^ h^jhiji^ 

{i,j},k i,j,k i j,k 

j 

Substituting this in (12321) yields (l2J0]) and (|2TT]) . □ 

Remark 2.7. Note that (|2.10p is equivalent to (j2.ip . Similarly, the first 
equation in Lemma 12.51 savs that = fjj, which is geometrically obvious. 
This provides an alternative, more geometric approach to Lemmas 12.61 and 
[231 

Lemma 2.8. The second variations of the face areas are given by 

A-i — ^ ^ hijhijif. 

Proof. This follows from Lemma 12.51 since hij = hijk = because of their 
linearity in h. □ 

Lemma 2.9. The second variation of the volume is given by any of the 
following two formulas. 

(2.13) YoV =^hiAi. 



(2.14) 



i 
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Proof. By differentiating (j2.10p and taking into account hi = 0, we ob- 
tain (f2l3]) . Differentiating (ETTDi yields 



Vol" =Iyi + ^ ^ 



2 ^ 2 

Combining this with (j2J3]) yields (pJij) . □ 

Remark 2.10. Another way to write the first and the second variations of 
Vol is 

Vol = - ^ hihijhijk = - ^ hihijhijk = - ^ hihijhijk, 

i,j,k i,j,k i,j,k 

Vol — ^ ^ hihijhijf^ — ^ ^ hihijhijf; — ^ ^ hih^jhiji^. 

i,j,k i,j,k i,j,k 

2.5. Proof of Theorem 12.11 The proof of Theorem 12.11 is based on the 
following key lemma. 

Lemma 2.11. Let h be such that Ai = for some i. Then we have 

(2.15) Ai < 0. 

Besides, equality in (I2.15P holds only if £ij = for all edges bounding the 
face Fi. 

This lemma will be proved later in this subsection. 

Proof of Theorem \2.1\ We prove Theorem 12.11 in a reformulation given in 
Subsection 12. II if h is such that = for all i, then hi = (aji'i) for some 
a G M^. 

If = for all i, then (12131) imphes 

(2.16) Vol" = ^/liii = 0. 

i 

On the other hand, by Lemma [2.11l we have Ai < 0. Besides, without loss of 
generality we can assume that hi > for all i (just choose the origin € M'^ 
inside Q). Hence (I2.14p implies 

Vol" = ^hiAi < 0. 

i 

By comparing this with (I2.16p . we deduce that Ai = for all i. Then, by 
the second part of Lemma [2.111 we have iij = for all edges of Q. It follows 
easily that h translates the polyhedron as a rigid body, thus hi = {a, Vi) for 
some a € M". Theorem 12.11 is proved. □ 

Fix an index % and consider the function : [/ ^ R. It can be written as 
a function of (hij). Here j varies over all faces of Q adjacent to the i-th face. 
Without loss of generality assume j = 1 , . . . , m and introduce new variables 

9j = hij, j = l,...,m. 
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Formally speaking, we consider a linear map 
(2.17) ^ R"", 

(hi) H> (gj) = (hij). 
The map factors through $ and induces a map 

A:V^R, 

where V = ^{U) C R™ is a neighborhood of g^ = ^{h^). We have 

m 

i=i 

where ij = iij are linear functions of g. 

By our construction, gj = hij are the support parameters of the face Fj 
with respect to the projection of to span(Fi). We have 

ogj ogjdgk dgk 

Lemma 2.12. The second differential D^A of the area function on the space 
of convex m-gons with fixed edge directions has signature (+i, O2, — m-s)- 
Besides, D^A at a point g takes a positive value on the vector g. 

Proof. First, let us show that 

(2.19) dimkeri:»2y4 = 2. 

By (j2.18p . D^A coincides with the Jacobian of the map g ^ I. Thus 

(2.20) gGkerL>2^ ^=> t = Q. 

But £ = implies that g is induced by a parallel translation of the polygon. 
Such deformations form a 2-dimensional space, and (j2.19p follows. 

Consider the space Qm of all convex m-gons. A convex polygon is deter- 
mined by its edge normals (/ij) and support parameters {gj)- The numbers 
{gj) must satisfy a system of linear inequalities (with coefficients depending 
on (/ij)) expressing the fact that all edge lengths are non- negative. For every 
collection (/Zj) of m different unit vectors positively spanning M^, the collec- 
tion {gj = 1) satisfies this system, as (/i,(7) corresponds to a circumscribed 
polygon. Thus Qm retracts to the space of configurations {^j) of m unit 
vectors positively spcinning M . It follows that Qm, 

is connected. 

Due to (|2.19|) . the rank of D^A is constant over the space Qm- As Qm 
is connected, it follows that the signature of D^A is also constant. Thus in 
order to determine the signature at the point we need, it suffices to compute 
it at a point we like. 

Let us compute the matrix of D^A- Let ajj^i,j = l,...,m, be the 
exterior angle between j-th. and {j + l)-st side {j + 1 taken modulo m). 
Then the formula on Figure [3] implies 
di ■ 

(2.21) — ^ = —(cot aj-i,j + cot Ojj^i), 
ogj 



(2.22) 
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Note that the matrix of D^A does not depend on g. This is because A is a 
homogeneous polynomial of second degree in (gj). By the same reason we 
have 

D^Aig,g) = 2A>0, 

which proves the second statement in the lemma. 

To determine the signature of D^A, put ctjj+i = Then we have 



D'A 



sin^ 

m 



(—2 cos 
1 




2ZL 

m 



-2cosi? 




1 

-2cos^ 



\ 



The spectrum of this matrix is 



2(cos ^-cos 21) 



k = l. 



-2 cost/ 



m > , which con- 



tains exactly one positive (for k = m) and two zero [k = 1 and k = m — 1) 
eigenvalues. The lemma is proved. □ 

Proof of Lemma \2.11\ Change the variables for the function Ai with the 
help of the map ()2.17p . Let g = d<&(/i) be the variation of g induced by h. 
By assumption, we have 

j 

On the other hand, since (ij) are linear functions of g, we have 

^dgk 



which results in 



D^A{g,g) = Y,Qf^9jgk 



^gjij 



3 



0. 



That is, the vectors g,g & are mutually orthogonal with respect to the 
symmetric bilinear form D^A. By Lemma 12.121 g is a positive vector for 
D^A, and D^A is negative semidefinite on the orthogonal complement to g. 
Thus we have 

A = D^A{g,g) <0. 

If i = 0, then g G keiD'^A, which by I^^^Wil) implies Ij = for all j. The 
lemma is proved. □ 

Remark 2.13. In general, the second differential does not behave well under 
change of variables. However, it does in our case, because the map (|2.17p is 
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linear. Therefore the second variation Ai in Lemma 12.111 is associated with 
the second differential D^A from Lemma 12.121 

Ai = D''A{d^h),d^h)) = D^A{g,g), 

which was implicitly used in the above proof. 

2.6. The case of a non-simple polyhedron Q. In the course of our proof 
we assumed that the polyhedron Q is simple, that is each of its vertices be- 
longs to exactly three faces. This assumption implies that the combinatorics 
of Q is preserved when its support parameters vary slightly. 

If the polyhedron Q = Q{h^) is non-simple, then Q{h) may be combinato- 
rially different from Q even for h close to /i". A neighborhood C/ of /i" S 
is subdivided into cells (U^) with A indexing simple combinatorial types 
of perturbed polyhedra. On each cell, the function Vol is a third degree 
homogeneous polynomial in (hi): 

\ol{h) = V^{h) for h G U^, 

so that function Vol is piecewise polynomial in a neighborhood of h^. It can 
be shown that 

Vol G C^{U), Ai E C\U). 

But Ai may fail to be C^, which is bad because second differentials of face 
areas play a key role in Subsection 12.51 

This problem can be resolved as follows. Choose any simple combinatorics 
A from a neigborhood of h^, and work with and Af' instead of Vol 
and Ai. Geometrically this means that we view our non-simple polyhedron 
Q{hP) as a member of a family of simple polyhedra, with some edge lengths 
vanishing; by varying in U, we allow negative edge lengths to appear. 
One can check that all arguments in Subsection 12.51 go through. 

This modification proves Theorem 12. II for non-simple polyhedra. 

3. Metric rigidity of convex polyhedra 

3.1. The theorem. Let P C M'^ be a convex polyhedron with vertices 
pi, . . . ,Pn- Assume that P is simplicial, that is all of its faces are triangles. 
Subsection 13. 81 explains how our arguments change in the non-simplicial case. 

Definition 3.1. An infinitesimal deformation of P is a collection of vectors 
Qi E M^, i = 1, . . . ,n. Each qi is thought of as a vector applied at the point pi. 
An infinitesimal deformation (qi) of a polyhedron P is called isometric, if 

(3.1) \\{Pi+tqi)-{pj+tqj)\\=0, 
"■'^ t=0 

for all edges piPj of P. In other words, if lengths of all edges don't change 
in the first order as vertices (pi) move linearly with velocities (qi). 

A simple computation shows that condition (j3.ip is equivalent to 

{pi -Pj,qi - qj) = 0. 
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Definition 3.2. An infinitesimal deformation {qi) of P is called trivial, if 
the map pi i— )• qi is the restriction of an infinitesimal isometry of M'^ ; 

= Api + h, Ae so(3), b G M^. 

An infinitesimal isometry of ffi^ preserves in the first order distances be- 
tween all pairs of points; therefore every trivial infinitesimal deformation of 
P is isometric. 

Definition 3.3. A polyhedron P is called infinitesimally rigid, if every iso- 
metric infinitesimal deformation of P is trivial. 

Theorem 3.4 (Legendre-Cauchy-Dehn). Every convex polyhedron in M'^ is 
infinitesimally rigid. 

This is usually referred to as Dehn's theorem. Let us explain why we 
prefer a different attribution. Cauchy [10] proved a global rigidity statement: 
two convex polyhedra with the same combinatorics and pairwise isometric 
faces are congruent. His proof was based on ideas presented by Legendre 
in \31\ note XII, pp. 321-334]. Dehn was the first to state and prove the 
infinitesimal rigidity theorem in [15]. However, in the footnote on the first 
page of [59], Weyl remarks that the argument in [10] carries over word by 
word to yield a proof of Theorem 13.41 and reproaches Dehn for not citing 
Cauchy. Note also that the Cauchy's "arm lemma" is more immediate in 
the infinitesimal context. 

3.2. A reformulation. Without loss of generality we may assume that the 
coordinate origin € M'^ lies in the interior of P. For every vertex pi of P, 
put 

r?=\\p.\\- 

The length of an edge pipj will be denoted by 

^ij = \\pi -PjW- 

Triangles OpiPj, where piPj ranges over all edges of P, cut the polyhedron P 
into triangular pyramids. These pyramids have a common apex at 0, their 
bases are faces of P. Let r € M" be a point close to = (rf)^^^. Change 
lateral edge lengths of the pyramids from to while keeping base edge 
lengths £ij fixed. A metric space glued from the new collection of pyramids 
(by the old gluing rules) is called a warped polyhedron. A warped polyhedron 
is not embeddable in in general, because the sum coi of all dihedral angles 
at an edge Opi may be different from 27r, for some i. Denote by 

Ki = 2tT — IjJi 

the curvature of the warped polyhedron at the edge Opj. This yields a C°°- 
map 

U M", 
r I— )• At, 

where U C is a sufficiently small neighborhood of r^ . 
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Lemma 3.5. A convex polyhedron P is infinitesimally rigid if and only if 

' dKi 



dim ker 



drj 



3. 



dt 



Proof. We will establish a correspondence between isometric infinitesimal 
deformations of P and elements of the kernel of J!^ = ■ 

Let q be an infinitesimal isometric deformation of P and let p\= Pi + tqi 
be the corresponding linear motions of the vertices. Put 

(3.2) r\=\\p\l t\^ = \\p\-p)\\. 

This defines a family of warped polyhedra, this time with non-constant 
metric on the bomidary. By Definition 13. II and since the polyhedron remains 
embedded in M^, we have 

- 4 = 0, K* = 
i=0 

for all i, j. Modify the deformation (j3.2p by putting l\j = Although the 
functions might be not identically zero anymore, their time derivatives at 
t = will not change: 

^ij — ^ij =^ = 0. 

It follows that r E ker J^. 

In particular, if q is an infinitesimal rotation around an axis through the 
origin, then we have rj = ri, and thus r = 0. If g is a parallel translation, 
then r 7^ is a non-trivial element in ker . One can show that translations 
in linearly independent directions in produce linearly independent vectors 
r G R". Hence trivial infinitesimal deformations of P give rise to a 3- 
dimensional subspace of ker . 

In the opposite direction, let us associate to r S ker an isometric 
infinitesimal deformation of P. Let P1P2P3 be a face of P. Choose a vector 
Qi collinear with Opi so that when pi moves with the velocity qi, its distance 
from changes with the speed ri: 

. Pi 

Choose q2 collinear with the plane 0piP2 so that to satisfy the conditions on 
variations of r2 and ^12: 

g2,-fr\)=r2, {pi-p2,qi-q2) =0- 



And choose qs so that to satisfy the conditions on variations of ra, £13, 
and £23: 

(3-3) ( q3, 7p-{7) = h, {Pi - P3,qi - qa) = 0, {P2 - P3,q2 - q3) = 0. 



Now proceed to an adjacent face and determine q^ by conditions 

similar to (|3.3p . Continue the face path (pipiPi^i) around the vertex pi 
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until it closes at piPkP2- On the vector there will be four conditions, the 
fourth being the stability of ik2- It is easy to see that this condition follows 
from the first three due to 



so that is well-defined. In a similar way we can find qi for all i, all closing 
conditions being satisfied due to k = 0. 

Note that we made some voluntary choices at the beginning, with qi 
and q2, but starting from q^ everything was forced. It follows that the 
infinitesimal deformation q is determined by r uniquely up to an infinitesimal 
rotation. 

As we observed, trivial infinitesimal deformations of P generate a 3- 
dimensional subspace of ker J^. If there is a non-trivial isometric infini- 
tesimal deformation, then the corresponding r € ker lies outside this 
subspace, as the map q ^ r identifies only deformations that differ by an 
infinitesimal rotation. Thus dim ker = 3 implies infinitesimal rigidity. 
Vice versa, if dim ker > 3, then there is an r G ker not coming from 
a trivial infinitesimal deformation. Thus it determines a non-trivial iso- 
metric infinitesimal deformation q. That is, infinitesimal rigidity implies 



3.3. The discrete Hilbert-Einstein functional. In Section [3.21 we intro- 
duced warped polyhedra that have curvatures (kj) along their radial edges 
Opj. Recall that the boundary edge lengths {(.ij) are kept constant, while 
the radial edge lengths r = (rj)^]^ are allowed to vary in a neighborhood U 
of a point G M". Denote by P{r) the warped polyhedron with radial edge 
lengths r. Let Ajj be the exterior dihedral angle at the boundary edge PiPj. 
Note that {Xij), as well as (kj) are functions of r. 

Definition 3.6. The discrete Hilbert-Einstein functional of a warped poly- 
hedron P{r) is 




dim ker 



= 3. 



□ 



(3.4) 




By X^jj j} we denote the sum over all unordered pairs of i and j, so that 
every boundary edge is taken once. 



Lemma 3.7. We have 



(3.5) 



5 HE 



dri 



Lemma 13.71 can be derived from (and is equivalent to) the Schlafli for- 
mula for Euclidean polyhedra. We will prove it by a different method in 
Subsection 13.51 
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Corollary 3.8. A convex polyhedron P is infinite simally rigid if and only 
if the second differential of the Hilhert- Einstein functional has corank 3 at 
the point r"; 

dimker(L>^HE)|^=rO = 3. 

Proof In view of Lemma 13.71 the matrix of the second differential of HE 
is the Jacobi matrix of the map r i— )■ k. Thus Corollary 13.81 is simply a 
reformulation of Lemma 13. 5[ □ 

3.4. The total curvature of a warped spherical polygon. In a warped 
polyhedron P{r), consider the spherical link Si of the vertex pi. By defini- 
tion, Si is a complex of spherical triangles obtained from trihedral angles 
at Pi by intersecting them with a unit sphere centered at pi. In our case, 
these triangles are glued cyclically around a common vertex, forming a total 
angle of 27r — . The angles at the boundary vertices of Si are equal to the 
dihedral angles of P{r), that is to vr — Xij. The lengths of radial edges in Si 
are equal to pij , where pij is the angle at the vertex pi in the triangle Opipj . 
See Figure HI Denote 

Sij — COS Pij . 




Figure 4. The triangle OpiPj and the spherical link Si of 
the vertex pi. 

A complex of spherical triangles similar to Si (a set of triangles glued 
cyclically around a common vertex) will be called a warped spherical polygon. 
Note that Si is in addition convex, that is the angles at its boundary vertices 
are less or equal tt. 

Definition 3.9. The total curvature Ki of a warped spherical polygon Si 
on Figure\^ right, is defined as 

Ki — l^i ^ ^ Sij Xij . 

j 
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Lemma 3.10. We have 

(3.6) RE = ^riK,. 

i 

Proof. The orthogonal projection of to the hne piPj sphts the edge piPj in 
two segments of lengths r i cos pij and rj cos pji. Thus we have 

(3.7) £ij — f^i^ij ~t~ TjSji. 

(If the projection of lies outside the edge, the formula remains valid.) By 
substituting (13. 7p in the definition of HE, we obtain 

HE = ^ riKi + ^ {nsij + rjSji)Xij = ^ r^Ki + ^ ViSijXij 

i {i,j} i i,j 



and the lemma is proved. □ 

3.5. Variations of HE and Ki. Let r ^ U, and let r G be an arbitrary 
variation of r. Denote by 



HE- = A 

dt 



HE(r + tr) 

t=o 



the derivative of the function HE in the direction r. Similar notations Ki, 
Sij, Xij are used for the directional derivatives of other functions of r. 

Lemma 3.11. The first variation of the total curvature of the spherical link 
Si of the vertex pi in P{r) is given by 

Ki = Sij Xij. 
j 

Proof. The warped polygon Si is made of spherical triangles, see Figure HI 
Consider a triangle with edge lengths pij, pik- Its third side has a constant 
length, as it equals to the angle PjPiPk on the boundary of the polyhedron 
P{r). By applying Lemma lA. 21 we obtain 

(3.8) Oijk + PijSij + 'jikSik = 0, 

where Oijk is the angle between the pij and pi^ sides, and /3jj, 7jfc are the 
two other angles. As we have 

aijk = 271 - Ki, Pij + 7jj = vr - Xij, 

jk 

summing (13. Sp over all triangles yields 

(3.9) + ^ii^ij = 0- 
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It follows that 

j j 

and the lemma is proved. □ 
Lemma 3.12. For every edge PiPj, we have 

Proof. This follows from ([32]) and iij =0. □ 

Lemma 3.13. The first variation of the Hilbert- Einstein functional of a 
warped polyhedron P{r) is given by 

(3.10) HE' = ^^ 



Proof. Follows by differentiating p.6p and applying Lemmas 13.111 and 13.121 

HE' = i-jKi + nkj = fi{Ki + sjjXjj) sijXij 



— ^ ^ ''i'^i ~l~ ^ ^ ij'iSij -\~ ^jSji ~\~ r^Sij + — ^ ^ T'il^i- 

□ 

Lemma 13.71 follows from Lemma I3.13( as (13. 5p is just a reformulation 
of (13TI)]1 . 

Recall that the second variation of HE is the quadratic form in (f-j) asso- 
ciated to the second differential of HE: 

HE =o2he(,,,) = ^ _ 



HE(r + tr). 



Lemma 3.14. The second variation o/HE is given by any of the following 
two formulas. 

(3.11) RE" =^riki. 

i 

(3.12) RE- =Y,r^Y,^iA^r 

i j 

Proof. Equation ()3.1ip is a direct consequence of (j3.10p . In order to prove 
(|3.12p . transform (|3.1ip by using (|3.9p and Lemma [3.121 

HE = y ^ riki = — ^ij^ij ~ ~ {'f'iSij + '''j^jij^ij 
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□ 

3.6. Proof of Theorem 13.41 The proof is based on the following lemma. 
Lemma 3.15. Let r be such that ki = for some i. Then we have 
(3.13) 

j 

Besides, equality in (I3.13|) holds only if Xij = for all j . 

The proof of this lemma requires a detailed study of the second variation 
of Ki with respect to the variables (sjj). This is done in Subsection l3. 71 which 
ends with the proof of Lemma 13.151 Now we prove Theorem 13.41 assuming 
the validity of Lemma 13.151 

Proof of Theorem \3.4\ We prove theorem 13.41 in the reformulation obtained 
in Lemma [3.51 It suffices to show that if r is such that /ij = for all i, then 
r is induced by a trivial isometric deformation of P. 
If = for all then (|3.1ip implies 

HE" = ^riki = 0. 

i 

On the other hand, equation (I3.12P and Lemma 13.151 imply 
(3.14) HE" =^ri^s,jAij <0, 

i j 

because > for all i. Thus we must have 

SjjAjj = 

j 

for all i. By Lemma 13.151 this happens only if Ajj = for all edges PiPj. 
This means that the isometric infinitesimal deformation that corresponds to 
r does not change the dihedral anges of the polyhedron P in the first order. 
It follows that this infinitesimal deformation is trivial. □ 

3.7. The second variation of the total curvature of a warped spher- 
ical polygon. In Subsection [331 we defined warped spherical polygons and 
their total curvature. Warped spherical polygons appeared as spherical links 
Si of vertices of a warped polyhedron P{r), see Figured Consequently, the 
total curvature of Si was viewed as a function of r E M". 

In this subsection, we study warped spherical polygons on their own. We 
preserve the notations on the right of Figure HI but suppress the index i. 
Without loss of generality, assume that the index j numbering the radial 
edges varies from 1 to m. As before, lengths of boundary edges are assumed 
constant, so that a warped spherical polygon is determined by m parameters 



Sj = cos pj, j = l,...,m. 
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with s varying in a neigborhood V oi E M™. The point has the property 

k{s^) = 0, 

where k (a former Kj) is the angular defect of the interior vertex of the 
warped polygon. 

Consider the total curvature of a warped spherical polygon as a function 
of s: 

K:V M, 

(3.15) K{s) = K + ^SjAj. 

j 

Lemma 3.16. The first and second variations of K as a function of s are 
given by 



(3.16) i^ = 5]s,A,-, 



J 



(3.17) K = J2^Ar 

j 



Proof. Equation (j3.16p follows from Lemma 13.111 as the form of the first 
variation does not depend on the choice of variables. (Actually, in the proof 
of Lemma [3.111 we deal with coordinates (sij), so we implicitly used the 
invariance of the first variation.) 

Equation (j3.17p is a direct consequence of (|3.16p . Compare this with 
(I31T]1 . □ 

We have abused notation: in ()3.17p . double dots denote the second vari- 
ation in the variables {sj), while before we used them to denote the second 
variation in variables (r^). Don't ever try to substitute ()3.17p in (I3.12p ! 

Lemma 3.17. The second differential D^qK of the function (I3.15P at the 

point s^ has signature (+i, O2, — m-s)- Besides, the associated quadratic form 
is positive on the vector s"; 

(3.18) D%K{s\s^)>Q. 
Proof. First let us show that 

(3.19) dimkerL»2(o)i^ = 2. 
Lemma 13.161 implies 

— - A hence - ^ 

dsj dsjdsk 

so that the matrix of D^K coincides with the Jacobi matrix of the map 
s X. Thus 

s G ker D'^K A = 0, 
that is if and only if the variation s of radial edge lengths induces a zero 
variation of angles between boundary edges. It is not hard to show that at 
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s = such variations form a 2-dimensional space, namely they come from 
moving the interior vertex while the boundary of the polygon remains fixed. 

Let Vm be the space of all convex spherical m-gons with a marked interior 
point (i. e. convex warped spherical m-gons with «; = 0). It is easy to see 
that Vm is connected. Every element of Vm can be viewed as a point in a 
space of warped polygons with a fixed boundary metric, thus has a matrix 
D^oK associated with it. By the previous paragraph, the rank of D^oK is 
constant over Vm- Since Vm is connected, and D^^K depends continuously 
on an element of Vm, the signature is constant as well. Thus in order to 
compute the signature of D^qK for all warped spherical m-gons with k = 0, 
it suffices to do this for one such m-gon. 

Let us compute the matrix of D^qK. Let be the angle between two 

consecutive radial edges. Then we have by Lemma lA.ll 



cot a 



+ cot a 



sin^ Pj 



dsj 



dX 



1 



sm Ojj+i sm Pj sm pj+i 

2 f^^ „ . . . _ 27r 

s' 



We compute the signature of D^qK for ajj+i 
In this case 

/-2cos^ 1 

1 



and Pj 



DlK 



1 



sm — sm p 

m. • 







-2 cos 
1 



2n 





1 



-2 cos 



2n 



p G (0, f 

1 \ 





-2 cos 



2n 



This matrix has spectrum 



2 (cos ^ 



sin — sin^ p 



,k— 1, . . . 



m 



which contains 



exactly one positive eigenvalue (for k = m). The first part of Lemma 13.171 
concerning the signature, is proved. 

The positivity of the form D^qK on the vector follows from Lemma 
EH] below. ^ □ 

Let C §^ be a convex spherical polygon with a distinguished point po 
in its interior. The polar dual S* of S is the intersection of all hemispheres 
centered at vertices of S. We define the Euclidean polar dual of S as the 

at pq. See 
then the 



projection of S* from the center of to the plane tangent to 
Figure m Note that if some of the distances pj are bigger than ^, 
point po lies outside the spherical and Euclidean polar duals. 



Lemma 3.18. Let S be a convex spherical polygon with a marked interior 
point Pq. Let (s^) be distances from pQ to the vertices of S. Consider the 
space of warped spherical polygons obtained from S by varying (sj) while 
keeping the boundary edge lengths fixed. Then we have 



DioK{s^,s") = 2ATea{S^) 



26 RIGIDITY OF SURFACES AND THE HILBERT-EINSTEIN FUNCTIONAL I 




Figure 5. Fragments of S, S*, and S^. 

where 5g is the Euclidean polar dual to S, see the definition before the 
lemma. 

Proof. Draw the perpendiculars from po to the sides of S*. They cut the 
spherical polygon S* into quadrilaterals, each with a pair of opposite right 
angles (quadrilaterals may be self- intersecting) . The signed lengths of the 
perpendiculars equal f — /Oj, with angles ajj+i between them. The Euclidean 
polygon has a similar decomposition, with the same angles ajj+i but 
with perpendiculars of lengths cot pj, see Figure O Formula on Figure [3] 
gives us the lengths of the other two sides of a + l)-quadrilateral, so 
that we can compute its area: 

1 / cot Pj+i — cot Pj cos ajj+i cot pj — cot Pj+i cos ajj+i 
— cot Pj : '■ \- cot pj-^^l 



2 \ sinajj+i sinajj+i 

By summing over j and performing simple transformations, we obtain 

Area(5|) = ^Jsl = ^I),^oi^(.°, s% 

and the lemma is proved. □ 

Proof of Lemma \3.15[ Suppress the index i and consider k = ni and \j = Xij 
as functions of {sj) = (sij). Due to invariance of the first differential, upper 
dots on the left hand side of (j3.13p can be viewed as variations with respect 
to the variables {sj), and we find ourselves in the setting of the present 
subsection. 

By assumption, we have k = for the variation s at the point ■s'^. Together 
with equation (13. 9p (or (I3.16P ) this implies 



J2s^^Xj = -k = 0. 



On the other hand, 



Ah, J 



j j,k j,k 

Thus we have 

DlK{s^,s) = 



RIGIDITY OF SURFACES AND THE HILBERT-EINSTEIN FUNCTIONAL I 27 



which means that the vectors s and s are mutuahy orthogonal with respect 
to the symmetric biUnear form D^oK. By Lemma [3.17| is a positive vector 
for D^qK, and D^qK is negative semidefinite on the orthogonal complement 
to s^. It follows that 

DloK{s,s) < 0. 

Since, by (|3T7D . 

j 

the inequality ()3.13p in Lemma [3.151 follows. Equality ^jSjXj = means 

that s € ker D^aK. This implies Xj = for all j, as established in the first 
lines of the proof of Lemma 13.171 

The lemma is proved. □ 

3.8. The case of a non-simplicial polyhedron P. First of all, if P 
has some non-triangular faces, our definition of an isometric infinitesimal 
deformation (the second half of Definition 13. ip has to be modified. If we only 
require stability of lengths of edges, then, say, the cube would be considered 
as infinitesimally flexible. There are two possibilities. Either one requires 
that each face of P is moved by (qi) as a rigid plate. Or, one subdivides 
each face into triangles by non-crossing diagonals and requires stability of 
lengths also for the diagonals. The latter class of infinitesimal deformations 
is a priori larger than the former, and can be shown to be independent of 
the choice of subdividing diagonals. 

We choose the second possibility. Thus, P = P(r^) can be viewed as 
a simplicial polyhedron with some dihedral angles equal to vr. As we vary 
r in a neighborhood U of r", these dihedral angles can become less than 
vr. Undaunted by this, we carry out our arguments for non-convex warped 
polyhedra as well. As we come to warped spherical polygons, we consider 
also non-strictly convex and "slightly non-convex" ones. Everything goes 
through, including Lemma [3.171 For the positivity of the form D^o^^ on the 
vector we need only positivity of the area of the corresponding Euclidean 
polar dual. This is fulfilled, as the polygon S is a convex (albeit non-strictly) 
spherical polygon. 

Thus a simple modification of the definitions allows to extend the argu- 
ment to non-simplicial polyhedra. 

4. Duality 

4.1. Duality between second derivatives of Vol and HE. There is an 
apparent duality between constructions and arguments in Sections [2] and [31 
This is the same kind of duality as between Cauchy's proof of Theorem 13.41 
in [To] and Alexandrov's proof of Theorem 12.11 in [U Chapter XI] . Cauchy 
studies a deformation of the spherical link Si of a vertex pi, marks with -|- 
or — the j-th vertex of Si if the angle (Xij) increases or decreases during 
the deformation and shows that either at least four sign changes occur as 
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one goes along the boundary of Si, or no vertex of Si is marked, i. e. the 
dihedral angles at pi are stable. Then a combinatorial argument involving 
the Euler characteristic of a polyhedron shows that there are no markings 
at all. Alexandrov considers variations of edge lengths along the boundary 
of a face Fi, marks lengthening and shortening edges, and follows Cauchy's 
argument. 

In this paper, new similarities show up. Deformations of Q in Section[2]are 
governed by support parameters (hi), while deformations of P in Section [3] 
are governed by radii (rj). There are functions Vol(/i) and HE(r) with similar 
variational properties, and infinitesimal rigidity is proved by showing that 

dimkerp^ Vol) = 3 and dimker(L»2 HE) = 3. 

The following lemma strengthens this analogy by making it quantitative 
rather than qualitative. 

Lemma 4.1. Let P C M.^ be a convex polyhedron with vertices pi, . . . ,pn 
and such that lies in the interior of P. Let Q = P* be the polar dual of P, 
that is 

Q = {x eM.^ \ {x, Ui) <h"i,i = l,..., n}, 

where 

Pi Pi ^0 _ 1 



— — , h- — 
\\Pi\\ * r?" 

Then we have 

(4.1) Vol = HE . 

Proof. A direct computation [6l Section 3] shows that 

92 HE _ cot ^ij + cot jij _ Vol 
dridrj iij sin pij sin pji dhidhj ' 

52 HE -r-^ cot j3ii + cot 7i,- Vol 

orf ^ ' lij sm pij sm pji dhf 

Here angles Pij and "fij were defined in the proof of Lemma 13. IH ipij is the 
angle at in the triangle QpiPj , and (.ij and pij are as on Figure [H The 
lemma is proved. □ 

An immediate consequence of Lemma 14.11 is that each of Theorems 12.11 
and 13.41 implies the other one (e. g. one can replace Subsections 13.31 to 13.71 
by Lemma l4.ip . On the other hand, a proof of identity (14. ip by a direct 
computation is more a question than an answer. There is no obvious reason 
why (j4.ip should hold, for example the functions HE(r'^ + a;) and \o\{h^ +x) 
are by no means equal. 

The next subsection partially "demystifies " identity (14. ip by pointing out 
a close relationship between HE and Vol in the context of polar duality in 
the 3-sphere. 
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4.2. Duality in spherical geometry. Let P C be a convex spherical 
polyhedron (i. e. the intersection of a finite number of hemispheres that 
contains no pair of antipodal points). Define the polar dual of P as 

P* = {xe8^ \ {x, y) <0 for ah y G P}. 
(We put §3 = {2; G m4 I ii^ii ^ g^j^^ ^.^ ^Yie scalar product in W^.) 
Then P* is also a finite intersection of hemispheres (take one hemisphere for 
every vertex of P). We assume that dimP = 3; then P* contains no pair 
of antipodal points. Let piPj be an edge of P. In P* there is a dual edge 
lying in the intersection of hyperspheres polar to pi and pj. Denote by iij 
the length of the edge PiPj, and by i*j the length of the dual edge. 

Theorem 4.2 (P. McMullen, [33]). For every S-dimensional spherical poly- 
hedron P and its polar dual P* the following identity holds: 

(4.2) Vol(P) + ^Y1 + Vol(P*) = vr^. 

The proof is based on Lemma 14.31 below. 
We have 

where Xij is the exterior dihedral angle of P at the edge PiPj, and A*^ is 
the exterior dihedral angle of P* at the dual edge. The first equation holds 
because the distance between pi and pj is equal to the exterior dihedral 
angle between the hyperspheres polar to pi and pj. The second holds since 
(p*)* = p and PiPj is dual of its dual. It follows that the summand in the 
middle of (|4.2|) is the discrete total mean curvature of dP and at the same 
time the discrete total mean curvature of dP*. 

More generally, let P be a convex polyhedron in the d-sphere S'^ for 
arbitrary d > 1, and let P be a face of P. We define the dual face of P* as 

F^ = {x eP* \ {x, y) =0 for aU y € F}. 

We have dimP-*- + dimF = d — 1. Every polyhedron is considered to be its 
own face and also to have as a (— l)-dimensional face. Then, clearly 

0-L = p*, = 0. 

Define the norm of a fc-dimensional face as its A: -volume divided by the 
/c-volume of the fc-dimensional sphere: 

As Volo(S'^) = 2, the norm of a vertex equals ^; we also put ||0|| = 1. 

Lemma 4.3. For every convex d- dimensional polyhedron P C S'^ the fol- 
lowing identities hold: 

(4.3) j;||F||-||F^|| = l, 

F 
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(4.4) • IIF^II =0. 

F 

Here the summation extends over all faces of P, including and P itself. 

Proof. Denote by F x F-^ the convex hull of F U F-^ in S'^. For every F, 
this is a ti-dimensional spherical polyhedron, and we have 

||F X F-^ll = ||F|| • \\F^t 

since span(F) and span(i<'^) are mutually orthogonal subspaces of M'^"'"^. It 
is not hard to see that the sphere decomposes as union of polyhedra 

= \J{F X F^), 

F 

with disjoint interiors. Thus we have 

Vold(§^) = J]Volrf(FxF^) 

F 

which implies equation (j4.3p . 

Equation (j4.4p is proved in a similar way, by replacing P* with —P*. The 
family of polyhedra F x — -F"*", each counted with multiplicity (— 1)'^™^, 
forms a covering with total multiplicity over every point of S"'. Therefore 

^(_l)dimFyQ^^(^ X F^) = 0, 
F 

and equation ()4.4p follows. Figure [6] illustrates the case d= 1. □ 




Figure 6. To the proof of Lemma [Ol 

Proof of Theorem^^ Since Yoh{^^) = 2ti'^ and Voli(Si) = 27r, for d = 3 
formulas (j43]l and (j44]l yield 

Vol(P) + ^ Area(5P) + + \ Area(5P*) + Vol(P*) = 2t^^ , 

Vol(P) - J Area(9P) + \Y^ ^i/ij - j Area(9P*) + Vol(P*) = 0. 

{hi} 
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By adding and subtracting these two formulas, we obtain 
Vol(P) + ^ E vr,- + Vol(P*) = vr^, 

(4.5) Area(9P) + Area(5P*) = Att 

Theorem 14.21 is proved. □ 

Remark 4.4. Equation (j4.5p is the discrete Gauss-Bonnet theorem for sur- 
faces in S^. Equation (j4.3p is the Steiner formula for the volume of the 
^-neighborhood of P. Equation ()4.2p is an analog of the Weyl tube for- 
mula. 

Lemma [4. 31 is proved by Peter McMullen in [33]. Section 4 of |33j contains 
also references to smooth analogs of formulas (j4.2p and (j4.5p . The main 
reference is Herglotz [22] (see also \2A\ Subsection 6.4]) who mentions that 
Fenchel has proved formulas 1 4 . 2 1 and (j4.5p by induction on dimension. There 
seems to be no written account of Fenchel's proof. 

Remark 4.5. Points of P* are poles of great spheres disjoint from the 
interior of P. This gives Theorem 14.21 an integral-geometric interpretation: 
a random great sphere intersects P with the probability 

^(vo1(P) + ^EM^'0- 

\ {iJ} / 

See also [43l Chapter 17, §5, Note 1]. 

Remark 4.6. Milnor on the last page of [35] suggests another way of prov- 
ing Equation (|4.2p . by deforming P to a point and integrating the Schlafli 
formula. It would be interesting to see whether the Schlafli formula can be 
derived from (|4.2p . 

4.3. Gauss and metric rigidity for convex spherical polyhedra. The 

following are the spherical analogs of Theorems 12.11 and 13. 4i 

Theorem 4.7. Let Q C S"^ be a simple convex polyhedron, and let {Qt) be 
a deformation of Q that keeps all dihedral angles constant in the first order. 
Then Qt is congruent to Q in the first order. 

This theorem extends to non-simple convex spherical polyhedra, if one 
requires stability of dihedral angles between the planes of any two faces 
having a vertex in common. A less restrictive generalization allows change of 
combinatorics and appearance of negative edge lengths, see Subsection 12.61 

To substantiate the analogy between Theorems 12.11 and 14.71 note that 
directions of face normals of a Euclidean polyhedron determine its dihedral 
angles, and that dihedral angles of a simple polyhedron determine angles of 
its faces, and thus, for a spherical polyhedron, face areas. 
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Theorem 4.8. Let P C be a simplicial convex spherical polyhedron, and 
let (Pt) be a deformation of P that keeps all edge lengths constant in the 
first order. Then Pt is congruent to P in the first order. 

A generalization can be stated for non-simplicial polyhedra by subdividing 
non-triangular faces, see Subsection 13.81 

Note that Theorems 14.71 and 14.81 are equivalent: just put Qt = P^ . 

One can prove Theorem 14.81 following the arguments used in Section [3] 
to prove its Euclidean analog. Warped spherical polyhedra are defined in 
an obvious way. In place of the Hilbert-Einstein functional consider the 
functional 

S{P) := 2Vol(P) + Y,riKi + 

« {hi} 

where P = P{r) is a warped polyhedron with radii (r^). Schlafli's formula 
implies 

dS _ 

~ — — Kj-i . 

dri 

Functional S can be described as a discrete gravity action with non-zero 
cosmological constant, cf. [Ml Subsection 6.5]. 

When we want to copy the approach of Section [2] to prove Theorem 14. 7^ 
it undergoes more substantial changes and becomes more close to the ap- 
proach to the dual theorem. Let qo be an interior point of Q. By dropping 
perpendiculars q^qi to the faces and then perpendiculars qiqij to the edges 
of (5, we cut Q into "bricks". A brick is a polyhedron combinatorially equiv- 
alent to the cube and with six right angles at the edges not incident to the 
vertices qo and qijk ■ A brick is uniquely determined by lengths hi , hj , of 
edges adjacent to q^ and by dihedral angles at the edges adjacent to qijk 
(just choose qo as a unique point at distances hi^hj^h^ from the planes of 
a given trihedral angle). We vary "support parameters" {hi) of a brick de- 
composition of Q while leaving dihedral angles of Q constant. As a result, 
singularities around the edges q^qi appear. Denote by t^i the curvatures at 
q^qi and put 

S*{Q) :=2Vol(Q) + ^/iiKi. 

i 

Then again Schlafli's formula implies 

dS* _ 
dhi 

As deformations of Q that keep dihedral angles constant can be identified 
with vectors h such that /i = 0, Theorem 14.71 reduces to a statement about 
the kernel of D'^S*. 

Lemma 4.9. Let P be a warped spherical polyhedron, and P* be its polar 
dual. Then we have 

S{P) + S*{P*) = 2tt^. 
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Proof. Let be a spherical cone manifold obtained by taking a warped 
product of [0, vr] with the spherical link 5*0 of the warped polyhedron P. 
(The manifold is a 3-sphere with 7i singular meridians of curvatures (^i); 
meridians are arranged as the cone singularities in the link of the point 
Po & P-) Both warped polyhedra P and P* can be embedded in S^. 
Similarly to Lemma 14.31 and Theorem 14.21 we can show that 

the right hand side being half the volume of S^. On the other hand, radii of 
P and heights of P* are related by /ij = vr — rj. Therefore 

S{P) + S*{P*) = 2 Vol(P) + iiji*^ + 2 Vol(P*) + + hi)Ki 

{iJ} i 

= 7r(27r — Kj) + vr Atj = 27r^, 

and the lemma is proved. □ 

Thus in the spherical case the duality between Gauss and metric rigidity of 
convex polyhedra is perfect, and their proofs that use the second derivatives 
of S* and S are simply the same. 

4.4. Duality in hyperbolic-de Sitter geometry. Consider the hyper- 
boloid model of the hyperbolic space 

= {x G I ||x||3,i = -l,xo > 0}, 

where ||x||3^i = —Xq + x\ + x\ is the Minkowski scalar product. The 
orthogonal complement x-^ C does not intersect H^, but it intersects the 
one-sheeted hyperboloid 

dS^ = {x G I ||x||3,i = 1} 

which is a model of the de Sitter space. For every convex polyhedron P C H"^, 
define its polar dual as 

P* = {xe dS^ U MiL I (x, y)3 J < 0}, 

where is the antipodal copy of H"^. In particular, dP* is a convex poly- 
hedral surface in dS^. Hyperbolic-de Sitter duality was used by Rivin [41] 
to study geometry of hyperbolic polyhedra. 

The Minkowski scalar product induces a semi-Riemannian metric on dS'^, 
and the boundary of P* is space-like with respect to this metric. Between 
lengths and angles in polar objects in H'' and dS'^ there is a correspondence 
similar to that in S"^, see Subsection 14.21 As a consequence. Gauss (respec- 
tively, metric) rigidity of a convex polyhedron in 'M? is equivalent to metric 
(respectively, Gauss) rigidity of its polar dual in dS'^. 



34 RIGIDITY OF SURFACES AND THE HILBERT-EINSTEIN FUNCTIONAL I 



A variational approach to the metric rigidity of convex hyperbohc poly- 
hedra is based on the discrete gravity action 

SiP) :=-2Vol(P) + ^ 

where P is a warped hyperbohc polyhedron. Similarly to the previous Sub- 
section, we have 

S{P) + S*{P*) = 0, 

where S* is the gravity action without the boundary term. This is implied 
by the following analog of Theorem 14.21 

Theorem 4.10. For every convex polyhedron P C and its polar dual 
P* C dS'^ U Hi the following identity holds: 

(4.6) Vol(P) -^Yl ^^/u- + Vol(P*) = 0. 

Note that P* is the union of a whole hyperbolic space Mi and of an 
infinite end of d§^. However, there is a consistent way to define a finite 
measure Vol(P*), [ails]. 

Similarly to Milnor's approach in the spherical case. Theorem 14.101 is 
proved in jSOj by integrating the Schlafli formula. 

Remark 4.11. Formula (j4.6p has the following integral-geometric interpre- 
tation: the (motion-invariant and appropriately normalized) measure of the 
set of all planes that intersect a convex body P C equals 

{iJ} 

In particular, this quantity is monotone under inclusion and always positive. 

Integral non-euclidean geometry is dealt with in [431 Chapter 17]. A 
proof of Theorems 14.21 and 14.101 indicated in |431 Chapter 17, §5, Note 1] is 
by obtaining first a similar formula in the smooth case and then going to 
the limit e for the boundary of e-neighborhood of P. 

4.5. Shearing vs. bending. There is another kind of duality between The- 
orems 12.11 and 13.41 relating metric rigidity of P with Gauss rigidity of P, 
rather than with that of P* . 

An infinitesimal isometric deformation of P can be described by assigning 
to each face of P an infinitesimal isometry of M^, so that these isometrics 
agree on the edges. (Faces move as rigid plates joined by hinges.) An 
infinitesimal isometry of is an infinitesimal screw motion, that is a vector 
field 

^{x) = T] X X + T, ri,T e M^. 
Let = {i]i,Ti) be the infinitesimal isometry associated to a face Fj. Let 
Fij = Fj n Fj be an edge of P. The condition (,i\Fij = CjlPij is equivalent to 

iVi ~ Vj) X X = Tj — Ti for all x € Fij, 
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which imphes that the vector rji — rjj is parallel to Fij : 
(4.7) rji — rjj || Fij for all edges Fij. 

Now, instead of rotating each face Fi according to the vector r/j, let us 
translate the plane of Fi by r/j (consider a plane as a set of points rather 
than as a geometric figure, so that the horizontal component of r/j cannot 
be neglected). By (14. 7p . the two translations of span(i<ij) differ by a vector 
parallel to Fij. This implies that the boundary of Fi is translated together 
with Fi, except that some new edges may appear. Namely, if Fi and Fk had 
only a vertex in common, the translations can split this vertex and create 
an edge between Fi and F^. If there is no vertex splitting at all, then the 
polyhedron P is translated as a rigid body. 

Even if vertex splitting happens, the area of each face changes by o(t), 
if we do translations by tr]i. It follows that putting hi = {r]i,Ui) we obtain 

= for all i, in the notation of Subsection 12. 1[ 

Thus to every isometric infinitesimal deformation of P there corresponds 
an infinitesimal Gauss image preserving deformation of P. The correspon- 
dence can be inverted, by solving the above equations for (ij). 

This observation was made by Weyl in [59j who was inspired by a work 
of Blaschke [5] where a similar correspondence for smooth surfaces was in- 
dicated. See also [U Chapter XI, §3]. A similar correspondence is known for 
ideal hyperbolic polyhedra, where it is sometimes described as the duality 
between shearing and bending. Indeed, with (dj) viewed as rotations, the 
difference di — dj represents the change of the dihedral angle between Fi and 
Fj, whereas if [di) are viewed as translations, then di — dj is the shift of Fi 
along Fij relative to Fj. 

4.6. Statics and polarity. Infinitesimal rigidity of a polyhedron is equiv- 
alent to static rigidity of its 1-skeleton (with diagonals of non-triangular 
faces added) , see [42l [26] . Basically, statics deals with duals of vector spaces 
of isometric, respectively trivial, infinitesimal deformations. 

In [BOl Theorems 3.1 and 3.2], Whiteley establishes a correspondence be- 
tween statics of a 3-dimensional polyhedron and of its polar dual. The 
polarity can be taken with respect to any quadric. This reflects the projec- 
tive nature of static (and hence infinitesimal) rigidity, see next subsection. 

4.7. Projective invariance and infinitesimal Pogorelov maps. Infin- 
itesimal rigidity of a bar-and-joint framework in is invariant under pro- 
jective transformations. This fact follows from a projective formulation of 
statics, see [111 [26]. 

Projective invariance provides a link between infinitesimal rigidity in W^, 
S*^, and W^: an infinitesimally rigid framework in remains infinitesimally 
rigid when viewed as a framework in a projective model of S*^ or W^. The aris- 
ing correspondences between infinitesimal isometric deformations are called 
infinitesimal Pogorelov maps, see [38[ [26] . 
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5. Miscellaneous remarks 

5.1. Regge action and Steiner formula. The sum ^^TiKi appearing in 
the definition of the discrete Hilbert-Einstein functional for warped poly- 
hedra makes sense for an arbitrary compact closed manifold glued from 
Euclidean simplices. It was introduced in |40j and is sometimes called the 
Regge action. It was shown in [12] that rim converges to a constant 
times the total scalar curvature if a sequence of piecewise Euclidean mani- 
folds converges in some good sense to a smooth manifold. 

The sum 2 j} ^ij^ij discrete analog of the total mean curvature. 
It appears as a coefficient at in the formula for the volume of an e- 
neighborhood of a convex polyhedron (as the total mean curvature does for 
convex bodies with smooth boundary). This was noticed by Jakob Steiner 
[39]. The convergence of "^^i jj ^ij^ij to a constant times the total mean 
curvature follows from the continuity of mixed volumes with respect to the 
Hausdorff distance, see [38] . 

5.2. Mixed volumes and the Alexandrov-Fenchel inequality. Our 

constructions in Section [2] are closely related to the theory of mixed volumes, 
see e. g. [H]. In particular, the proof of Theorem 12.11 is copied from (a part 
of) the proof of Alexandrov-Fenchel inequalities, [H [3l UH]. The relation 
between mixed volumes and derivatives of the volume is apparent from the 
formula 

Vol(/i + tk) = Vol(/i) + 3t Vol(/i, h, k) + St^ Vol(/i, A;, k) + \o\{k) 

that may serve as the definition of mixed volumes Vol(/i, /i, k) and Vol(/i, k, k). 
(One has to assume that polyhedra Q{h) and Q{k) are combinatorially iso- 
morphic, in order that addition of support parameters h and k correspond 
to Minkowski addition of polyhedra.) Right hand sides of the formulas in 
Remark 12.101 give alternative expressions for mixed volumes. 

To prove Theorem 12. H it is sufficient to show that the symmetric bilinear 
form Vol has corank 3. The proof of Alexandrov-Fenchel inequalities 
goes further. It establishes that D^Vol has signature (+i, O3, — „,_4), and 
in general shows by induction that the signature is {+i,Od, —n-d-i) for d- 
dimensional polyhedra with n facets. 

5.3. Existence theorems. Infinitesimal rigidity can sometimes be used to 
prove existence theorems. Assume that an object is infinitesimally rigid with 
respect to some parameters. If the space of objects and the space of param- 
eters have equal dimension, then the inverse function theorem provides us 
with local existence and uniqueness of an object with given parameters. In 
other words, the map {object} {parameters} is a local homeomorphism. 
A local homeomorphism with good topological properties is a global home- 
omorphism, which implies that any given set of parameters defines a unique 
object. 

This method is used in [4] to prove the following two theorems. 
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Theorem 5.1 (Minkowski theorem). Let ui,V2, ■ ■ ■ ,i^n G IR^ be unit vectors 
spanning M.^, and let Ci, C2, . . . , Cn be positive numbers such that 

i 

Then there is a unique convex polyhedron in with outer face normals (vi) 
and respective face areas (Ci). 

Theorem 5.2 (Alexandrov theorem). Let g be a Euclidean cone-metric 
on ( every such metric can be obtained by gluing a set of Euclidean tri- 
angles). Assume that the angles at all cone points are less than 2tt. Then 
there is a unique convex polyhedron in such that g is an intrinsic metric 
of its boundary. 

Note that in both theorems the combinatorial structure (which pairs of 
faces are adjacent, which pairs of cone points are joined by edges) is not 
given in advance and is practically impossible to determine without finding 
the polyhedron in question. 

Properties of the functionals Vol and HE suggest a variational approach 
to both theorems. Indeed, 

^ (^o\{h) - hiC^ = Mh) - a 

implies that the polyhedron in the Minkowski theorem is a critical point 
of the function Vol(/i) — ^ • hiCi on the set of polyhedra with outer face 
normals {vi). Similarly, 

9HE _ 

implies that the polyhedron in the Alexandrov theorem is a critical point of 
HE on the set of warped polyhedra. As stated in Subsection 15. 2| the second 
differential Vol has only one positive eigenvalue. By a happy coincidence, 
Vol is concave on the hyperplane {^ihiCi = 1}. At the maximum point 
(whose existence can be shown by a simple trick) face areas Ai are propor- 
tional to Cj, so that a scaling yields the desired polyhedron. This proof is 
due to Minkowski himself, [36], see also [U Chapter VH, §2]. 

The situation is more complicated with HE and the Alexandrov theo- 
rem. The discrete Hilbert-Einstein functional on warped polyhedra is nei- 
ther concave nor convex, moreover the signature of its second differential is 
non-constant. Nevertheless, similarly to Subsection 13.71 one can show that 
kerD^HE = {0}, provided that Kj > and the areas of Euclidean po- 
lar duals of the vertex links are positive. This allows to apply the inverse 
function theorem to the map r 1— >■ k and to find the desired polyhedron by 
constructing a family of warped polyhedra with k ^> 0. This proof is given 
in [6]. 

Yuri Volkov, a student of Alexandrov, studied the discrete Hilbert-Einstein 
functional with the aim to find a variational proof of Alexandrov's theorem 
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similar to that of Minkowski theorem. His proof in [53\ I56j) doesn't use the 
function HE, but is similar in the spirit. The polyhedron is found by mini- 
mizing the sum of radii (r^) over all warped polyhedra (this time with apex 
at a boundary vertex) with negative (kj). See also [55] which is reprinted in 
the Appendix to where Volkov studies HE in order to give some a priori 
bounds for the embedding problem. 

In some theorems of Alexandrov or Minkowski type, the functional HE, 
respectively Vol, happens to be concave. This is the case with the Alexan- 
drov convex cap theorem [25] , and with the Alexandrov and Minkowski- type 
theorems for convex hyperbolic cusps [18^ I17j . However, in all these cases 
proofs of existence of a maximum point are quite difficult. 

5.4. Related work. Our proof of infinitesimal rigidity of convex polyhedra 
in Section [3] is related to works of Pogorelov and Volkov [39\ [5l] and of 
Schlenker [36]. In [461 Section 3] an alternative proof of Lemma [3.171 can be 
found. 

Duality between metric rigidity of P and Gauss rigidity of P* involving 
theory of mixed volumes appears in Paul Filliman's work [19]. Lee general- 
ized this duality [30] to higher dimensions, relating it with Peter McMullen's 
theory of weihghts on polytopes [34] . See also the work of Tay, White, and 
Whiteley [5l]. 

The approach of Section [3] was extended by Schlenker to star-shaped 
polyhedra with vertices in convex position, [17|- A further development of 
this technique is given in [28| . There the signature of HE for an arbitrary 
triangulation of a convex polyhedron is determined. 

5.5. Directions for future research. 

Problem 1. Can equation ()4.ip be obtained as a limiting case of Theo- 
rem |li2]? 

Problem 2. Give a proof of Theorem 14.101 similar to McMullen's proof of 
Theorem 14.21 Probably, all prerequisites are contained in [13]. See also |22| . 

Lemma 14.11 relates metric rigidity of P and Gauss rigidity of P* , Sub- 
section |4?6] relates metric rigidity of P and P* , and Subsection 14.51 relates 
Gauss and metric rigidity of P*. In the smooth case, this forms a part of 
Darboux's wreath, [241 Subsection 5.3]. Sauer j44j found a discrete analog 
of Darboux's wreath for polyhedral surfaces with quadrilateral faces and 
four-valent vertices. 

Problem 3. Is there a Darboux's wreath for polyhedral surfaces of arbitrary 
combinatorics? 

Problems HHT] are discussed in more detail in Subsections 11.3111.51 

Problem 4. Reprove the Calabi-Weil rigidity of compact closed hyperbolic 
3-manifolds by showing 



(5.1) 
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where n is the number of vertices of a triangulation of a hyperbohc manifold. 

Problem 5. In a similar way, prove the infinitesimal rigidity of compact 
hyperbolic 3-manifolds with convex polyhedral boundary. 

Problem 6. In a similar way, reprove the infinitesimal rigidity of compact 
hyperbolic cone-manifold with all cone angles less than 2tt. 

Problem 7. Reprove Cheeger's vanishing theorem with a discrete Bochner 
technique. 

The next problem generalizes problems H] and [6] to higher dimensions. 

Problem 8. Define discrete Einstein manifolds in dimensions greater than 
3 and prove their infinitesimal rigidity under suitable assumptions. 

To approach the above problem, one can start from the following con- 
sideration. Gauss infinitesimal rigidity holds for convex polyhedra in all 
dimensions. Convex polyhedra of higher dimensions are also infinitesimally 
rigid, in fact they are too rigid. An alternative extension of the infinitesimal 
rigidity theorem to higher dimensions could be infinitesimal rigidity in the 
class of discrete Einstein metrics with a given restriction to the boundary. 
To define discrete Einstein metrics for a star-like triangulation of P, con- 
sider infinitesimal deformations h = u ol P* that preserve volumes of facets 
in the first order. Putting r = u should define an infinitesimal Einstein 
deformation of P. 

Problem 9. Define infinitesimal Einstein deformations of warped convex 
polyhedra in dimension d > 3 and prove Einstein infinitesimal rigidity of 
convex polyhedra. 

Appendix A. Some trigonometry 

The first two lemmas concern a spherical triangle with side lengths a, 6, 
c, and values a, /3, 7 of the respective opposite angles. 

Lemma A.l. The partial derivatives of the angles as functions of side 
lengths are: 

da _ 1 _ 1 
da sin b sin 7 sin c sin /3 ' 

da cot 7 

db sin b 

Proof. Straightforward calculation using spherical cosine and sine laws. □ 

Lemma A. 2. Let the side lengths b and c change with velocities b and c, 
while the length a remain constant. Then the corresponding variations of 
the angles satisfy the equation 

(A.l) d -|- /3 cos c -|- 7 cos 6 = 0. 
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Figure 7. Variation of the area of a spherical triangle when 
one of the side lengths is preserved. 

Proof. This can be proved by tedious computations using Lemma lA.ll We 
give an alternative geometric argument. 

Denote by area the area of the triangle. As area = a + /3 + 7 — vr, equation 
(jA.ip is equivalent to 

area = /3(1 — cos c) +7(1 — cos b) . 

The latter is proved on Figure [7] (recall that the area of a spherical cap of 
radius c equals 27r(l — cose)). □ 
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